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STIMMAR Y 

Modern aircraft design trends have been toward 
aircraft with long, dense fuselages, smaller aero- 
dynamic surfaces, swept wings and higher wing loadings. 
These aircraft are required to fly at high altitudes. 
Flight data from such aircraft have indicated that the 
linearized dynamic equations do not always predict the 
dynamic behavior accurately. In addition, certain cross- 
coupling effects were developing between the lateral 
and longitudinal motion of the aircraft. 

The effect of the trend toward concentrating the 
mass in the fuselage was investigated by making an 
engineering approximation of the change in the moments 
of inertia caused by this change in mass distribution. 
Four configurations were considered, each one approach-= 
ing nearer the limit of a missile or body of revolution, 
Since the lower lift curve slopes and smaller wing areas 
of the swept-wing aircraft sometimes require flight at 
very high angles of attack, the effect of inclination 
of the principal longitudinal axis to the flight path 
was investigated. The effect of sweep was approximated 
by varying the effective dihedral parameter, os 0 
In addition the effect of varying Ong , the "weather- 
cock-stability™ parameter, was studied since long fuse- 
lages affect this parameter adversely. 

The investigation was made using a Reeves Electronic 
Analogue Computer (REAC) by setting up a mathematical 


model of a representative military-type airplane flying 


= = 


= ~ @ oa te 


— re & _t«—- << 


ee ieee SS ie iin “el 


Sc oo om oe 


2 aaa — 2 ome! 
— a eS ee Se oe ah 

ec oe goa eemmecsime imal 

> | an ae Se 


7 





Zon 


at a Mach number of 0.6 at sea level and at 40,000 

feet. The mathematical model was based on equations 

of motion which were extended to include products of 

the perturbation velocity terms. The effect of changing 
the mass distribution, the inclination of the longi-~ 
tudinal principal axis, and variation of Ong and 

CL, was investigated by changing the mathematical 
model as required by the parameter change considered 

and obtaining the REAC response to a fixed or standard 
step-function control input. 

The results indicated that the lateral motion 
coupled into the longitudinal motion and that this 
effect was increased as the fuselage became longer 
and denser. Altitude aggravated this coupling as it 
did any trend toward instability caused by a change in 
the other parameters considered. The rudder caused 
greater coupling than the aileron. A positive inclina- 
tion of the principal longitudinal axis caused an 
increase in dynamic stability and a decrease in coupling 
while a negative inclination had the opposite effect. 

A decrease in Cn, produced a marked trend toward in- 
stability and an increase in coupling while a decrease 
in Ch had the opposite effect. 

The investigation required by this thesis was 
carried out in the Computer Center of the University 


of Minnesota, 
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INTRODUCTION 

Fighter aircraft of recent design show a trend 
toward long fuselages where most of the mass is con- 
centrated, smaller aerodynamic surfaces, high rates of 
roll, and increased use of sweepback. Flight data 
from such aircraft indicate that the linear equations 
of motion, which result after several simplifying 
assumptions, do not always predict accurately the 
dynamic behavior. In addition, these same design 
trends have required more emphasis to be placed on the 
consideration of the dynamic stability of an aircraft 
in the design stage. In earlier practice, dynamic 
stability «Stent took care of itself after the static 
stability requirements were met. 

The effect on coupling and dynamic stability 
caused by the parameters principally affected by these 
design trends was investigated by means of a mathe- 
matical model entered on the Reeves Electronic Analogue 
Computer (REAC). An attempt was made to use a representa- 
tive or typical aircraft of the fighter type. The 
mathematical model was based on equations of motion 
which were extended to include the products of the 
velocity change usually assumed to be zero. 

The mass distribution change was investigated by 
considering four configurations each more nearly 
approaching a missile from the standpoint of its moment 
of inertia properties. The effect of the inclination 


of the principal longitudinal axes was considered at 
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inclination angles of zero, plus ten, and minus ten 
degrees. The effect of Gn, and Coy. was studied by 
changing the mathematical model as required by a change 
in tnese parameters over a range of values. Altitude 
effects were considered by calculating the dimensional 
derivatives for 40,000 feet, changing the mathematical 
model accordingly, and running the program as it was 
done for the sea level case. 

The effect of a parameter change was observed by 
obtaining successive REAC responses to a standard step 
function control input. Only the parameter investigated 
was changed so that any variation in the time history 
response was attributed to the parameter change. The 
data was presented by superimposing these successive 
time histories. The trend of the effect of the para- 


meter change considered could then be easily observed. 
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DEVELOPMENT OF THE PROBLEM 

The classical approach to airplane dynamics is 
outlined here in a very brief manner. The airplane 
was considered to be a rigid body. Six equations 
were required to define its motion; one force equation 
and one moment equation for each axis. These equations 
included the force and moment components due to aero= 
dynamic and mass reactions together with the force of 
gravity and thrust effects. Second order effects were 
neglected and the equations simplified to a degree 
which permitted hand solution. 

Appendix I develops the equations of motion 
used in this investigation. The classical approach 
has been used in general but certain departures have 
been made. In Appendix II the equations developed in 
Appendix I have been simplified by assuming that the 
airplane initially is in straight and level, steady, 
gliding flight. The thrust effects were allowed to go 
to zero since they are not critically essjmtial ina 
dynamic stability study. 

Derivative terms of order higher than one in the 
series expansion of the expressions for the forces and 
moments were neglected. The effect of these terms is 
known to be secondary, data was available only for the 
first derivative, and inclusion of these terms, even if 
available, would require more components than the REAC's 


provided. A second assumption generally made is that 
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the product of the perturbation velocities is essenti- 
ally zero. This yields two independent sets of three 
linear equations. This assumption was not made. The 
equations simplified in Appendix II are presented below: 
(qa="7) a= Xyu + Xaq + Xue ur tXypur + VR-urg -g6 
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The axes used and the sign convention observed is 
presented in Table I and illustrated in Fig. 8 which 
heads Table I, 

Ref. 2 points out that motion of the longitudinal 
symmetric type does not induce motion of the lateral 
assymmetric type. An examination of the equation set 
presented above indicates that there are no terms which 
would couple symmetric motion into the assymmetric 
motion if the assymmetric motion were initially zero. 

Duncan in Ref. 2 also points out that when devia- 
tions are so small that their squares and products can 
be neglected, there will be no coupling of the lateral 


motion into the symmetric motion. When the deviations 
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are not small longitudinal motion will persist alone, 
but a purely lateral motion will, in general, induce 
longitudinal motion. The motion will couple into the 
forward velocity through the term vnr3 the pitching 
velocity through the terms %n (Ete), vn Ta and 
ile a * and into the normal velocity, through the 
term pvr .o It can be observed that the coupling terms 
are modified generally by moments and products of inertia. 
The study of the effect of the trend in mass dis- 
tribution toward a long, dense fuselage was made by 
considering four cases starting with a hypothetical 
fighter airplane of conventional type. The character- 
istics of this basic airplane are presented in Table II. 
The fuselage was considered to represent 60 per cent 
of the airplane's mass. It was taken to have a maximum 
width of five feet and to have an eight-to-one slender- 
ness ratio. The pitching moment of inertia calculated 
using these assumptions and the approximation that, for 
purposes of making an engineering estimate, the fuse- 
lage can be considered a long, slender cylinder resulted 
in a calculated value of the pitching moment of inertia 
which was very near the Ly of the airplane. As the 
mass became completely concentrated in the fuselage it 
was assumed that the slenderness ratio would increase 
to twelve to one, the fuselage width would remain the 
same, and that the pitching and yawing moments of 
inertia could be calculated using the long cylinder 


approximation. In this final form, the aircraft would 
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be a missile with ibe, equal to lig and i. very small. 
Figs 10 illustrates the engineering approximation 
that was made of the variation in moments of inertia 
as the fuselage became longer and denser. The mass 
distribution changes were made while the aerodynamic 
damping terms, control power terms, mass and other 
effects were considered constant. This would, of 
course, be unrealistic but it has the advantage of 
Studying the effect of a change in one parameter in 
the "partial derivative" sense. 
The product of inertia term resulted largely from 
the fact that stability axes were used. When the 
% -axis is so oriented in an airplane in the steady 
flight condition that it is parallel to the relative 
wind, the Eulerian axes are referred to as stability 
axes. NACA TR 589 indicated that the Il,. terms were 
unimportant for conventional aircraft, however, Ref. 7 
showed that the I,, term may have a pronounced effect 
on the lateral stability of some high speed airplanes. 
It has already been noted that the rs term is involved 
in the coupling of the lateral and longitudinal motion. 
If the longitudinal axes are inclined at an angle 
Y. to the flight path, the moments of inertia about the 
flight path axis and the axis normal to the flight 
path are: 
Te = Taco + Lp, aint 
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The angle ¥ largely results from the wing incidence 
angle and the angle of attack. 

The trend toward swept wings and long fuselages 
makes G, m4 One increasingly important. Long fuse- 
lages adversely affect Cu. ,» the "weather-cock" 
parameter. Swept wings have very large values of 
effective dihedral thus affecting, Ce, , the 
effective dihedral parameter, The effect of a change 
in these parameters on coupling and on dynamic stability 
was investigated. The change in the static derivative 
characteristics was selected for study since general 
flying characteristics are influenced primarily by the 
static stability parameters according to Ref. 6. 

The coefficients such as lL, are dimensi onal 
stability derivatives. They are direct functions of 
airframe geometry, mass, mass distribution and flight 
conditions; in addition to being direct functions of 
the non-dimensional stability derivatives such as, 

a » These basic non-dimensional stability derivatives 
are, in turn, implicit functions of airframe geometry, 
mass, mass distribution and flight conditions. The 

values of the non-dimensional stability derivatives 

for the basic airplane appear as Table II. The dimen- 
Sional stability derivatives are defined in Tables III 

and IV. The values of the dimensional derivatives were 
calculated and appear as Table V. Refs. 1, 2, 3, 4, 7, 

8, and 9 discuss these derivatives, their characteristics, 


and how they are determined. They are assumed to be 
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constant for a certain range of applicability and flight 
conditions. 

The dimensional stability derivatives can be used 
directly as numerical coefficients in the dynamic 
equations of the airframe based on real time. Thus, 
stability derivatives in this form are used in setting 
up the mathematical model for the analogue computer, 
The use of this form of derivative permits the exten- 
Sion of the non-dimensional derivatives to other flight 
conditions with only secondary error. This was done 
to obtain the mathematical model of the airplane at 
40,000 feet. Certain values of the non-dimensional 
derivatives which were specifically applicable to 
flight at this altitude were used while others were 
extended by the change in the values of _ and U , 

Phillips in Ref. 4 considered the effect of 
steady rolling on longitudinal and directional stabil~ 
ity. The equations which he used can be arrived at 
from the equations developed in Appendix I by applying 
the assumptions that he made in order to arrive at an 
equation set which could be handled analytically. 

He retained the product terms which couple the longi- 
tudinal and lateral motions. He found that rapid 
rates of roll caused instability if the directional 
and longitudinal stabilities are different when the 
rolling frequency exceeds the lower of the pitching 
and yawing natural frequencies of the non-rolling 


airplane. The instability existed only while the air- 
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plane was rolling. In airplanes of short span and 
high fuselage density flying at high altitudes, 
this instability may cause dangerous attitude 
changes. Further, he found that when an airplane 
rolls about an axis which is not alligned with its 
longitudinal principal axis, forces were introduced 
which tended to swing the fuselage out of line with 
the flight path. These are the inertial forces 
usually neglected when the linearized equations are 
used. 

In Ref. 4, Phillips made many assumptions in 
an effort to obtain equations which were capable of 
hand solution yet would predict the effect on coupling 
of steady rolling. In the investigation conducted 
incident to this thesis, the effects of the three motions, 
roll, pitch, and yaw, on coupling as produced by step-= 
function control inputs were investigated. As many 
influencing factors were retained as the REAC's compon- 
ents would permit. 

Solutions of the equations of motion used in this 
thesis would require the simultaneous solution of six 
differential equations, three of second order, made 
non-linear by product terms such as Gk and pr .e 
The labor involved wuld have been prohibitive. The 
equations were fitted to two REAC's and four servo- 
multipliers. The solutions were recorded simultane- 


ously on two four-channel Sanborn recorders. 
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EQUIPMENT AND PROCEDURE 

The equipment used to conduct the investigation 
consisted of two REAC's, four servo-multiplier units, 
two four channel Sanborn recorders, and a one-shot 
multi-vibrator time switch. The entire computing 
System is illustrated in the photograph presented as 
Fig. 2. The Sanborn recorder is shown in Fig. 3, the 
lateral motion and control function computer in Fig. 4, 
the longitudinal motion computer in Fig. 5, and the 
time switch in Fig. 6. 

Each of the four servo-mechanisms drove three 
multiplying potentiometers. These multipliers provided 
=> of the product~of the inputs. This restriction 
made it necessary to build up the perturbation velocities 
to ten times their actual value in order to obtain a 
Simple product suchas pr. The additional amplifi- 
cation required by the servo-multipliers made computer 
balance very delicate and this was a source of difficulty. 

The Sanborn recorder made time history records of 
four quantities simultaneously. The imprint was made 
by a heated pen on sensitized paper. The point did not 
touch the paper, thus eliminating the drag usually 
associated with a recording pen. 

The time switch was designed to provide an accurately 
timed step-function as a standard input. Basically, the 
switch is a one-shot multi-vibrator which fires just 


once when the operating button is pressed. This allows 
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current to flow in the coil of a relay in the plate 
circuit of the "normally-off" tube. The relay contacts 
send plus or minus 100 volts through the aileron, rudder, 
or elevator function switch, whichever is closed. 

The time that the relay remains closed is governed by 
the RC time constant controlled by the variable potentio-~ 
meter. A range of from approximately 0.5 to 2.5 

seconds was available. The electrical schematic is 
presented as Fig. 7. A photograph of the "bread-board" 
model appears as Fig. 6. The contacts of the relay 
carrying plus or minus 100 volts would appear electri- 
cally just ahead of the function switches in the REAC 
solution schematic which appears as Fig. 9. 

The six equations of motion were put in computer 
form in Appendix II. These equations are II-13 through 
II-18, A computer schematic for the simultaneous 
solution of these equation was laid-out and appears 
as Fig. 9. From the schematic, patchboards were wired 
for the lateral motion and control computer and the 
longitudinal motion computer. This division of the 
problem was based on the wiring scheme which would use 
the fewest amplifiers, The apparent symmetry was 
incidental. The numerical values of the quantities 
entering the problem were translated into potentiometer 


settings. The values recorded were U, 100, ww, 
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and loo on one recorder; Vy, 10d 10V; and Sen the 
other recorder. The plus and minus 100 volts was 
placed on one side of the timer relay contacts. The 
other side of the relay contacts was connected to each 
of the control function switches, plus or minus, which. 
ever was required by the schematic. 

The graph representing the estimated mass dis- 
tribution changes as the fuselage becomes longer and 
more dense is presented in Fig. 10. Values of the 
moments of inertia of three configurations were picked 
off this graph. New stability derivatives were cal- 
culated as necessary for these configurations. All of 
the moment derivatives such as L, were affected. These 
calculated values were translated into potentiometer 
settings and appear as Tables VII and VIII. 

The configurations investigated were called 
Case I, II, III and IV. The lateral computer was 
referred to as the "blue" computer and the longitudinal 
computer as the "red" computer. These are occasionally 
referred to in the tables using the abbreviations B and 
Re 

The computers were balanced and the recorders 
calibrated. The patchboards were then installed on 
the computers. A rudder, elevator, or aileron control 
movement was then selected and the recorder tape set in 


motion at the desired speed. Pushing the time switch 
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button sent an accurately timed standard step-function 
through the control switch selected. The solutions of 
the equations of motion in response to a step-function 
input were recorded as time histories of the variables 
on eight recorder channels simultaneously. The next 
case was then entered on the potentiometers and the 
solutions obtained as before, These solutions are 
presented as Figs. 220, 128, 134, 133, lipA, and 1B. 

The dimensional stability derivatives were recal- 
culated using data for 40,000 feet and the required angle 
of attack increase obtained from Table II. Where 
specific non-dimensional derivatives were not available 
for this altitude, those for low altitude were extended 
by using the values of Pp and VU at 40,000 feet in the 
formulae of Tables III and IV. These stability deriva- 
tives and other factors of the problem which changed at 
40,000 feet were translated into potentiometer settings 
and run off as before. The settings appear in Tables IX 
and X. The results for Cases I and IV are presented in 
Figs. 15A, 15B, 16A, 16B, 17A and 17B. It should be 
pointed out that the airplane which was marginal in 
dynamic directional stability at sea-level was unstable 
at 40,000 feet. In order to "fly" even marginally at 
this altitude the N, potentiometer was set to 0.0115 
and the ee potentiometer to zero. This had the effect 
of increasing the "“weather-cock" stability parameter 


about six times and reducing the effective dihedral 
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parameter to zero. The computer procedure was as 
before. It should be noted that solutions were 
obtained only for Cases I and IV modified as des- 
cribed above. 

The effect of inclination of the longitudinal 
principal axis to the flight path was investigated 
by calculating values of moments and products of 
inertis for inclination angles Y of plus ten, zero, 
and minus ten degrees. The equations below obtained 


from Ref. 6 were used: 
i = lly. Cos y sa Th. sin’ | 


ioe; i. Cos*y + Ty. sin? M 
Typ * ~(Ts. — i han ¥} Cos | 


It was found that only the 1. value varied significantly. 
New values of the coefficients containing [,, were 
calculated and translated into potentiometer settings. 
These are presented in Table XI. Several sign changes 
were required. These were accomplished by using inverters 
as necessary. The computer was operated as before 

using Case I at low altitude. 

The Ch and ct, programs were accomplished by 
changing these values over the range considered, 
recalculating the affected stability derivatives, and 
translating them into potentiometer settings. This 
was done for Case I at sea level and at 40,000 feet. 
These values are presented in Tables XII and XIII. 


The computer operation was similar to the other programs. 
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RESULTS AND DISCUSSION 


Time histories of the variables u,ur, 46,7, 

Pp , 6 , o » and y were recorded as the airplane 
responded to a step function control input for each 
change in a parameter being investigated. To study 

the trend of the effect of the parameter changes, the 
curves were superimposed, each curve being drawn with 

a characteristic line. These composite solution curves 
are presented as Figs. 12A through 23B. The standard 
step function time plot appears as Fig. ll. It should 
be noted that the A and B parts of the figures each 
represent one of the four channels which were simultan- 
eously recorded. Where separate lines for each para- 
meter change can not be observed, it should be assumed 
that the curve lies behind the previous curve, or curves, 
and that no significant change was caused by the new 
value of the parameter. 

The mass distribution program is presented in the 
Figs. 12A, 12B, 13A, 13B, 14A, and 14B. The elevator 
response curves show no significant changes caused by 
the changed moments of inertia. The elevator set-up 
the phugoid motion which remained essentially constant. 
The short semetiod motion was virtually dead-beat damped. 
The amplitude of Q in Case IV is reduced since the air= 
craft was more sluggish to pitch change due to the in- 
creased pitching moment of inertia. A slight dutch- 
roll was induced. This was not anticipated since the 


equations indicated that there would be no coupling 
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from the symmetric motion into the lateral motion if 
there was no initial lateral motion. A possible 
explanation might be that there was some unbalance 

in the servo~-multipliers or jitter in the lateral 
motion computer due to small line current variations. 
The same effect might be observed in flight tests 
caused by slightly unsymmetrical control application 
or small trim unbalance. Spiral instability was set- 
up presumably by the same mechanism which caused the 
Slight dutch-roll to start. 

The cross-coupling set-up by the aileron and 
rudder motion can be observed in thew plot. The 
coupling caused by the more assymmetric motion induced 
by the rudder step is several times that of the aileron 
cross-coupling. The aileron’s cross-coupling decreased 
as the values of I, and 1, were increased. This was 
unexpected but can possibly be explained by the fact 
that two opposing effects were taking place. The 
coupling increased but the airplane became more slug- 
gish since the pitch and yaw moments of inertia were 
increased. This sluggishness can be observed in the 
progressive decrease in amplitude and frequency of 
and V . The large decrease in roll moment of inertia 
had very little effect on ¢$ and b : 

The altitude effects on coupling and dynamic 
stability are presented in Figs. 15A, 15B, 16A, 16B, 


17A, and 17B. Many of the observations made in the 
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sea-level program apply here. The short period 
oscillation damping decreased increasing the period 

of this motion. The aileron coupling was larger for 
Case IV, that is, as the fuselage became long and dense. 
The same effect was observed in the coupling caused by 
the rudder. This time there was little amplitude 
difference between the aileron and rudder coupling. 
The sluggishness caused by increased We enidwl =. 
together with decreased aerodynamic damping, made the 
swings progressively larger in amplitude and period. 
The roll angle did not significantly change. 

In the program which considered the effect of 
changing the inclination of the longitudinal principal 
axis to the flight path, the 1. program essentially, 
only the rudder and aileron step input response was 
investigated. AS Yj became more positive, dynamic 
stability increased; more negative, dynamic stability 
decreased. At an Yi of minus ten degrees, dynamic in-~ 
Stability developed. From the coupling standpoint, the 
same effect was observed. There is less coupling when 
(| becomes more positive; more when \ becomes more 
negative. The rudder curves show these effects more 
clearly than the aileron curves. 

The effects caused by varying the "weather-cock" 
parameter, Cry 3 is presented in Figs. 20A, 20B, 21A 
and 21B, Only the aileron and rudder forcing functions 


were considered. The curves for two values of Cre 
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at high altitude were superimposed on the sea-level 
curves. At both altitudes similar effects were 
observed. Coupling increased as Cn, decreased. The 
dynamic lateral stability decreased as Cus decreased. 
This was to be expected, of course, since decreased 
tail arm, or tail area, lead toward lateral instability. 
As the angular velocity changes increased with this 
instability the coupling increased. 

The curves for the igi program are presented 
in Figs. 22A, 22B, 23A and 23B. The entire sea-level 
program was presented but only the curve for a 4, 
of 0.20 at 40,000 feet was added. As “a was increased 
dynamic, lateral instability quickly developed along 
with increasing coupling. The effect was more apparent 
in the motion caused by the rudder forcing function. 
The reduced damping at altitude hastened the develop- 
ment of instability and the coupling increased. The 
curve for a Shy OF 0.20 was marginally stable at sea 
level but unstable at 40,000 feet. 

The parameters C,,, and Cy, had very significant 
effects on dynamic stability. The effective dihedral 
parameter, Cis was very critical. This would appear 
to make analysis of the dynamics of the airframe very 
important in the design stage of swept-wing and delta- 
wing aircraft. The REAC provides a powerful tool for 


such an analysis. 
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The use of a standard step-function input 
provided a good means for studying the effect of a 
parameter change. However, this method was inadequate 
from the standpoint of bringing out any instabilities 
caused by cross-coupling of the type investigated in 
Ref. he 

It should be pointed out that the dynamic be- 
havior observed was that of one airplane, or variation 
of one airplane. Although effort was made to make 
the mathematical model investigated "typical™ , it 
would not be reasonable to suggest that all airplanes 
behave in the manner of the airplane investigated. 

On the other hand, the procedures used in the analysis 
were Straight~-forward and it is felt that the dynamic 
behavior observed might be considered reasonably 


typical. 
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CONCLUSIONS AND RECOMMENDATIONS 


Mass distribution changes had no significant 
effect on the motion of the aircraft in the 
longitudinal plane, 

The motion of the aircraft in the symmetrical 
plane coupled into the lateral motion in the 
form of a slight dutch roll and also triggered 
the spiral instability mode. However, this is 
believed to have been due to slight assymmetry 
in the computing system. 

The assymmetrical motion of the airplane coupled 
into the motion in the symmetrical plane and is 
most significant as a change in w . This would 
be equivalent to a change in the angle of attack. 
The coupling effect increased as the fuselage 
became longer and more dense. 

The assymmetrical motion caused by the rudder 
had a greater coupling effect than the motion 
caused by the aileron. 

At high altitudes both the coupling and the 
trend toward instability increased. 

At a positive inclination of the principal 
longitudinal axis to the flight path, there 

was less coupling and increased dynamic stability. 
A negative inclination angle had the opposite 


effect. 
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A decrease in the "weather-cock" stability 
parameter, Ch, , resulted in a decrease in 

dynamic stability and an increase in coupling. 

An increase in the effective dihedral parameter, 
Cr, , resulted in a decrease in dynamic 
stability and an increase in coupling. 

The standard step-function input which supplied 
the forcing function provided a means for studying 
the effect of a parameter change on the coupling 
and the dynamic stability of the airplane. How- 
ever, it does not appear adaptable to the investi = 
gation of instabilities which may develop from 

the coupling caused by rapid rolling or yawing. 

An investigation might do this by using a pro- 
grammed input which would reverse before the REAC 


limits were reached, 


The dynamic behavior observed in this investigation 


was that of the mathematical model of one airplane or 


variation of one airplane. Effort was made to use a 


mathematical model which could be considered reason- 


ably "typical", However, it would not be reasonable 


to suggest that all airplanes would respond similarly. 
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SYMBOLS 


initial airspeed, feet per second 

mass density of air, slugs per cubic foot 

wing span, feet 

wind area, square feet 

weight of airplane, pounds 

mass, slugs 

acceleration of gravity, feet per second per second 
aspect ratio 

Mach number 

wing chord tip, feet 


wing chord root, feet 


moment of inertia about principal longitudinal 
axis, slug-feet* 

moment of inertia about principal tranverse 

axis, slug-feet* 

moment of inertia about principal normal axis, 
slug-feet* 

moment of inertia about flight path axis, slug-feet* 

moment of inertia about the pitching axis, 
slug-feet@ 

moment of inertia about the yawing axis, slug-feet* 


product of inertia with respect to the flight path 


axis arid axis normal to flight path, slug~feet* 
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Tim 


force, vounds 


moment of momentum, (slug-ft)-ft per second 
Lift force 
q & 


rolling moment coefficient , 


lagt cochirci 6m , 


Rolling mement 
gq Sb 


lateral-force coefficient , tetera} force 

aerodynamic pitching moment coefficient about the c.g. 

drag coefficient , Drag force 

angle of bank, radians 

angle of pitch, radians 

angle of yaw, radians 

control deflection, radians 

angle of sideslip, radians 

flight path angle; angle between relative wind 
and the horizontal, radians 

angle between the longitudinal axis of the fuse~ 
lage and the relative wind, radians 

angle between the thrust line and the longitudinal 
axis of the fuselage, radians 

perpendicular distance from cg.to the thrust 
line, feet 

fraction of critical damping 

frequency, radians/sec. 

initial forward velocity in direction of positive 

axis (feet per second} 


initial side velocity (feet per second) 


initial normal velocity (feet per second) 
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forward velocity change after disturbance (feet 
per second) 

Side velocity change after disturbance (feet per 
second) 

normal velocity change after disturbance (feet 
per second) 

initial rolling angular velocity (radians per 
second) 

initial pitching angular velocity (radians per 
second) 

initial yawing angular velocity (radians per 
second) 

rolling angular velocity increment, radians per 
second 

pitching angular velocity increment, radians 


per second 


Aw yawing angular velocity increment, radians per 


Cp 


u 


second 

the change in drag coefficient with varying 
forward velocity 

the change in drag coefficient with varying angle 
of attack, per radian 

rate of change of lateral-force coefficient with 
rolling-angular-velocity factor, per radian 

rate of change of lateral-force coefficient with 
yawing-angular-velocity factor, per radian 

lateral force derivative, rate of change of lateral 
force coefficient with angle of sideslip, per 


racdien 
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the change in side force coefficient with variation 
in rudder deflection, per radian 

the change in side force coefficient with aileron 
deflection, per radian 

change in lift coefficient with variation in forward 
velocity, angle of attack and attitude remaining 
constant 

change in lift coefficient with varying pitching 
velocity, no change of airplane angle of attack 
as a whole, per radian 

the change in lift coefficient with varying angle 
of attack, "lift curve slope", per radian 

change in lift coefficient with changes in elevator 
deflection, per radian 

rate of change of rolling moment coefficient with 
yawing angular velocity factor, per radian 

effective-dihedral derivative, rate of change of 
rolling moment coefficient with angle of side- 
Slip, per radian 

damping in roll derivative, rate of change of rolling 
moment coefficient with rolling-angular-velocity 
factor, per radian 

change in rolling moment coefficient with variation 
in rudder deflection, per radian 

aileron effectiveness on aileron power, the change 
in rolling moment coefficient with change in 


aileron deflection, per radian 
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longitudinal static stability derivative, change 
in pitching moment coefficient with varying 
angle of attack, per radian 

the change in pitching moment coefficient with 
variation in rate of change of angle of attack, 
per radian 

pitch damping derivative, change in pitching moment 
coefficient with varying pitch velocity, per radian 

elevator effectiveness or elevator power, the change 
in pitching moment coefficient with change in 
elevator deflection, per radian 

directional stability derivative, weather-cock 
stability, rate of change of yawing-moment co= 


efficient with angle of sideslip (° 





an ee 
=) per radian 
damping in yaw derivative, rate of change of yawing 


moment with yawing angular velocity factor, 





: sC 
per radian ay 
2a 


rate of change of yawing moment coefficient with 
Nee 


rolling angular velocity factor, per radian Gee 

rudder effectiveness or rudder power, the change a 
in yawing moment coefficient with variation in 
rudder deflection, per radian 

change in yawing moment coefficient with change in 
aileron deflection, per radian 

change in fore and aft force with change in forward 
velocity, per second 


change in fore and aft force, change in pitching 


velocity, feet per second-radian 
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change in fore and aft force with change in 
normal velocity, per second 

change in fore and aft force with rate of change 
in normal velocity, per foot 

change in side force with change in yawing velocity, 
per radian 

change in side force with change in side velocity, 
per second-radian 

change in side force with change in rolling velocity, 
per radian 

change in side force with aileron deflection, 
per second~radian 

change in side force with rudder deflection, per 
second-radian 

change in normal force with varying forward velocity, 
per second 

change in normal force with varying pitching velocity, 
feet per second-radian 

change in normal force with varying vertical 
velocity, per second 

change in normal force with elevator deflection, 
feet per second*-radian 

change in rolling moment with varying yawing velocity, 
per second~-radian 

change in rolling moment with varying side velocity, 


per second-feet 
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change in rolling moment with varying rolling 
velocity, per second-radian 

change in rolling moment with rudder deflection, 
per second*-radian 

change in rolling moment with aileron deflection, 
per second“-radian 

change in pitching moment with varying forward 
velocity, per second-feet 

change in pitching moment with varying vertical 
velocity, per second-feet 

change in pitching moment with rate of change of 
vertical velocity, per foot 

change in pitching moment with varying pitching 
velocity, per second-radian 

change in pitching velocity with elevator deflection, 
per second*-radian 

change in yawing moment with varying side velocity, 
per second-feet 

change in yawing moment with varying yawing 
velocity, per second-radi an 

change in yawing moment with varying rolling 
velocity, per second-radian 

change in yawing moment with rudder deflection, 
per second*-radian 

change in yawing moment with aileron deflection, 
per second<-radian 

change in yawing moment with side-slip angle, 


per second*-radian 
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refers to initial steady flight condition 
pertains to the fore and aft axis 
pertains to the transverse axis 

pertains to the normal axis 


elevator 


aileron 


rudder 









Pie & 


Lp,¢ Longitudinal 
Or rolling OXI 


Lateral or 
pitching Axis Y) M,a,9 


Pe 5 
4 | Yaw ing K 
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Force ( parallel to 
axis) symbol x x 
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Ww 











Sign convention | Forward + 


Change in 


linear velocity (a 


Sign Convention Forward + Teweiant ce Downes 





| Moment about GhIS Rolling moment 
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Fermula M=+2pV*SeCy, N&© 2 eV'SbC, 

Sign convention rey 3) ae xy? re eae 
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Ral! Be - ; Pitch angle | Yow anne 
7 a 
Right roll + Nose uo + Nose right + 
Pitching velocity Yawing velocity 
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Right roth + Nose tup +. | Nose right + 
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LONGITUDINAL STABILITY DERIVATIVES 
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Dimensianar Derivatives or TRE Basic AIRPLANE 
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APPENDIX I 


Development of the Equations of Motion 


The equations of motion of the airframe are written 
by equating the forces and moments acting on the air- 
frame to the plane's reactions in accordance with 
Newton's laws. The reference axes are considered 
fixed in space. 

The aircraft is considered to be a rigid body; 
aeroelastic effects are not considered. The motion 
is referred to right-handed Cartesian axes. The axes 
used and the sign conventions observed are shown in 
Figs & and Table I. 

Newton's second law states that the rate of change 
of momentum of a body is proportional to the force 
applied to the body and the rate of change of the 
moment of momentum is proportional to the torque applied 
to the body. The mass of the airplane is considered 


constant for the duration of the dynamic maneuver 


considered. 

(1) 2 Fy = (rn) = om U 
> Fe a an V V 
2 fy a ( ) mm 
Fe = sb. (ren) * een We 

Zs) = any 

dat 

MA any 

at 

ZEN = dh 

at 
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Investigating the dynamics of an infinitesimal 
element of mass dm of the airplane referred to the 
Cartesian coordinate system yields equation set (3). 

(63) dl, « (aes) PAY, ~#xRdm - uy Qdm 

dhy - (+) Oum = xy Pam -ue Rdm 
Gag * (x*+ 4") Rdm - 2x Pam — Z24Qdm 

For a finite mass these become 

(4) hy s p | (t+ at)den : Q | x4 damn R(xedm 

hy Q (e+e) dm ~- R Jue am — Pfyxdm 
hes? Rl ue tyr) dma — P fex dm ~ Q fey dm 

The integral [cyt tet) dm ls recognized as the 

moment of inertia I of the whole airplane about the 


x-axis. Equation set (4) can be written 


h, Ply - "Gai “ae 
Ly + "Oly ~ Sane ey 
he® Rig Pile “i 
The equations of motion relative to inertial 


axes becomes: 


SF. may 

(5) ye mdi 
Lys m av 

7 at 

dF, « mM daw 

at 


(6) 2L= dime. = PI, 7s Pl, > Q Ty Boe, ~ ae -R lez 
2Mzé aby = ou, to Q yer R Lge erg, PL y, ~ Pile 
at e 8 e bd & a 
IN «dn, . 2 ie #RIZ- PL, Ply -@ ley ~QTey 
oe 
Equations (6) are written with respect to fixed 
inertial axes. It is expedient to use Eulerian axes 


with their origin at the c.g. which are fixed in the 


airplane and move with it. At any instant the airplane 





ibe 


in a maneuver has a motion referred to the Eulerian 
axeS and the axes have a motion referred to fixed 
space. The earth and its atmosphere are regarded as 
fixed space. The absolute acceleration of the aircraft 
referred to the Eulerian axes can be written 
(7} cee dV en 
= at 

where V_ is the total velocity. 

Similarly the absolute rate of change of the | 
moment of momentum is written 
(8) dhabs . dh en 


At at 
The first term arises from the change relative to 





the Eulerian axes fixed from instant to instant in the 
airplane and the second considers the change of the 
Eulerian axes relative to inertial axes fixed in space. 

Moments and products of inertia measured relative 
to Eulerian axes are independent of time since the mass 
is assumed constant for the duration of a maneuver. 
Many of the terms in equation set (6) go to zero for 
this reason. 

The term TH xV, of equation (7) can be written: 


—" 


(8a) wxVy = 


L 


$A PF! 


) 
Q 
V 


Cc ty z-b 


—_ —_ —_ 


Where taj, are unit vectors directed along the 
Eulerian X,Y, z axes. 

Expanding the cross product and combining with the 
components of the first term of (7) that equation can 


be written in component form as: 
ay = U +Qw -RV 
Az zw t Pyv- AU 
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By parallel reasoning (8) can be written: 
dh, 


(10) Gl, + dhe h@-h R= PL -@L4-RI4 +OR(I,-Iy) 
~ PQ Tra ~@ Tey RM yg, + PR Ty 





dh, dhy : 

a ie: +L. ee he Po. Q I, - RIu, ; Pile + PR( 1x-T,) 
-Qa Tey al ve p* Vi, ame ley 

dh 


at lab o GE + hy Pu . Q- R1e-Plye-Q1ay + PO (1y-1,) 
PR Dy leg HQ ene le 
With the Eulerian axes as chosen and shown in 
Fig. 8, the xz -plane is a plane of symmetry. There 
is both a positive and a negative value of y for 


each value of z ; consequently Tye + (42 dm=o and 


thee { xqam - a: 
The equations of motion of an aircraft referred to 
Eulerian axes can then be written: 
(Uy amma - om | U + Qw-rv] 
y Fy = m| V + RU- Pw | 
> Fe « mlLw +PV- Qu] 
PL* PT, = R ieee (a) rae 
Qty Ra - cae ee 
IN = Rig-Plve + PQ(14- Ly) +-ORTyve 


The left sides of the eqations in set (11) are 
the summations of the external forces and moments 


applied to the airplane in flight. These external 


forces are gravity forces, aerodynamic forces, and 
thrust forces. 
The gravity force may be considered to act at the 


center of gravity making no contribution to the 


summation of moments. 
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The resolution of the gravity force into components 
referred to disturbed Eulerian axes is shown in Fig. l. 

The force relations may be written 
(12) DF, = m[U +Qw- Rv] - Xs 

Fy + mV + RU- PW] - Y, 

DF, + am[W +Pv- Qu] - 2, 
by transferring the gravity terms to the right-hand 
side of the equation. The gravity terms X,, Y,;, and 
Z, are expanded in Fig. 1. 

The left side of all the force and moment equations 
now include aerodynamic and thrust forces as well as 
moments due to control surface deflection. 

The disturbed motion of an airplane at any instant 
can be considered the result of disturbing the airplane 
from some steady flight condition. This motion may 
be referred to as disturbed flight and the Eulerian 
axes as the disturbed axes. The forces acting on the 
airplane.under these conditions may be regarded as 
the sum of the force required to produce steady flight 
plus an increment required for the disturbance. Accord- 
ingly, each of the total instantaneous velocity com- 
ponents can be written as the sum of a constant velocity 
component during steady flight and a change caused by 
the disturbance: 

(13) Deter om 
V+ Ve + ar 


Wi = W, * uF 
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Substituting in (12) considering that es HB. om 
etc., and using the expansion of the gravity terms 
from Fig. 1 after setting terms like sinwsingso , 
sinb=Y and cos6+1 , the set of equations below 
can be written. The assumptions that the sine of an 
angle can be replaced by the angle, the cosine is 
equal to one and the product of sines can be replaced 
by zero produces results which, although approximate, 
are in good agreenent with physical results according 


to Ref. 7. These assumptions, of course, apply to the 


disturbed angles. 


(Uy) ZF tmLh + @QWet Weg +Qour - RoVe- RoV —Von turq~ TR 
+4 sin, - (g cw. sind.) p + (gcosO, cose) 6] 
DEL + mur + BVit Por +Vorp - ole - Quu-Ua-qut pr 
+ (goin Oo) 6 + (4 cobs 2 40) d (gem bs coh) | 
IF, ve ma LAr + US FR TUL + Rou - FoWr =F us Wee ~ up fee 
(Grind) y - 9 Ca Oo surb, — (genO, cng.) $ J 
De p Le-f lez t(QR, + Qn tRoq+qnr)(Iz-L,) 
- (PQ, + hq + Qept a) Tye 
2M: q L, +(PRo t+ Pon tRep +pnry(Ty - i) —(R+2eQntn) Tye 
+P. + 2Pop +r) Tye 
IN« Ale -p lez + (BQ tq +Qop rbq \(Ly-Lx) (QR Gon + Rogtng) Le 
The aerodynamic forces and moments acting on an 
airplane, the left side of the above equations, are 
functions of the flight conditions and of the deflection 


of the control surfaces. 
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Each of the forces (X\andz} and the moments 


(L,m,anaN ) can be expressed for small oscillations 
by expanding by MacLaurin's theorem neglecting second 


order terms to obtain: 


a 
Ke X, +u Rk ~ 4 oF + ur ee + Oe 


where 2x , Ls , etc. are partial derivatives of X 
eS 4 


with respect to u,q,etc. Exactly analogous 


expansions are made for Y,zZ,L,*1, and N . These 

forces and moments can be expressed: 

ag) %eX- Se eat wer + Ger. ss 
Y=y + ohn ts Bvt ety ss ee 


mee. Laer 354 Sh + SE 


Du 


. ig memes ach lene ag 2 


DR ov op Pe x5 = 
Mz Mo+iIMu+ 2 set ee pede ee <a 
R ul “e 2 Sy 
= 3230 + oN, + ON + ON 4 DN 
aa 
Or Dh sp P ae 5h OF 


Each of the terms in mie has a physical sig- 


nificance. X,,Y,andZ,,L,m. andN. are the forces 


aa 
and moments acting along and about the %*, 4, and z 
axes while the plane is in steady flight. The terms 
similar to 2X express the change in the given 
force or + caused by the given disturbance quan- 
tity. The coefficients of the form tL are the 
basic dimensional stability derivatives. A complete 


discussion of stability derivatives may be found in 


Rems. (17, (2), (37) (7), (8), and (9). 
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Thrust effects are part of the left sides of 
equation set (14). The equations for steady flight 


condition become: 


(16) Oe Tl cos & 
ee. = —T, Sita 
M, = Te 2; 


Since the Eulerian axes remain fixed with reference 
to the airplane during a disturbance, the thrust 


components relative to disturbed axes become: 





(17) X= T, cosé 
aati 
M = Ss 
where JT the thrust during the disturbance, is T, + AT « 
Al can be considered 
. rune > Fae tho 
and 
(19) Ko 1, case + Ceoss) BT y + (cess) 5 eos 
eer 
Z = ~ Te sun — (sin) ~ Cain) 2 2E eek 
Bape 


M= To; + 4 TK ae zs... 

The equations for steady flight are found by 
substituting the expressions obtained above for the 
thrust and the aerodynamic forces (and moments) into 
the left sides of equation set (14) and setting the 
disturbed quantities equal to zero: 

(20) X, - WsinQ. + To cose + O 
ae ae oe) 
7, + Wes. ~ Ls = 


exe 
Ua Ome ome. O 


Mo ot Ore. F O 
Not AO +t OS 


If now the complete equations for the disturbed 
airplane are written with both steady flight and 
disturbed terms they will be of the form: 


Aerodynamic forces + thrust forces = inertia force = 


gravity force 


Equation set (20) can then be subtracted from the 


complete equations of motion for the disturbed alr- 


plane mens: 
2X 2X 3X x Y 3T 
(ay) Bee tq 4 Shur + OR + eng )atu + (ms) Sem 


mi + Qo Mle + Weg + Cour turg = RM - Row -V, 
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>.) e 
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rp” CUTE 
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Dividing the force equations by m and the moment 


equations by the appropriate moment of inertia yields 


terms of the form: 


at. aa and | ®b 
Mm ou vary 
Replacing eM by Xu and { ay by L, simplifies 
xX 


the notation. These quantities are called dimensional 
stability derivatives. A list of these stability 
derivatives in terms of the dimensionless forms appears 
as Tables III and IV. The flight path angle 7, replaces 
QO.» The equations in final form can be written: 
(Ze) K (ex fae Kopiar toGaur + Coa § Meee + Gos ¥ Ts, Supe 
=u + Weg * Qour turg- Rov, ~Rnr ~ Vo -— wR 
— (9c i, nud.) v 
Ne + Yer Yep + V5, 5n1 YoqS, * F + Von Wap wp Mike An 
#Rout un - Pour -(g sino) p - g sim No sim 6, ~ (ger. corde] $ 
tou + 44 + £, our + Ee de “sing Tu ~ Taam tren ~ 
ir + PaVy + Pow + Nop - QU, ~Quu -Ve 4g ~Gu tp 
= (g Sur 6,) 6 + (G.crm@. an $,) d + (geo Gerd, = gG Coa, 
$n Tre + (QoW, +n + Req tgn)UTe=Ly)- (RO.t Pog, +Gup eq) ie 
= On +lyrtbhp-p + wre rls be : 
4 as ( Pol. + Pon +R-p +-pr) (1x-Iz) — (2 +2R ntnt) Txe 
+ (P3442 p+-p?) Txe = Ui dt +Mg4q + Mur wr Mal be 
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APPENDIX II 
Adaptation of the Equations of Motion 


to the Analogue Computer 


The dynamic equations of motion of the airplane 
were developed in Appendix I. These are summarized 


below: 
(a-1) Kyu eXyg + Xypur +X yur + cos§ are COM 


= i re Wo 4 - eae + Ur 4 ~ 1&Vo “Rove - Vo VR 
~ (gan to and,) a 
(m-2) Yah +Yv +Vpp + Ye 3, NE ar +U rn — Wop 
=wep + UR, UZ nn + Roy tun — Par 
(9 se) b — Gd, sndde — (gems cod) $ 
(ir-3) Aut 24.4 + Lup wr + &s 8,7 Sin gM =, Usp Bre © 
ut PRY, + Bw +Vop —QeoU, ~ Qu - oF Eu TP 
—(qon@o)O + (gre, br b,) > + (qa28, cer $,) ~ 9 ont, 
Ges) An n Tae +(Q, Qt Qer- oq + qr) ( — 4) 
Liza. + Poq. + Qop+ by) aXe 
Lin gon ee $+ Ls Om 


~y 


5) 4+ (Re Rn + Re pepal(Lacis) - (RE Hone) bee 
z 4 
+ (Po +2 Pop +p ) Uve = Ma ut Ma ay +M pwr + Me be 


(T-6) A-+ Tye +( PQ, +Peg+ Oop t Pq) (T Ty- Ix) 
<< 
+ Q, ‘ bent Beg +4qn) = Nap + Nt t+ Ne Ng Lie “ye 


These equations are general and allow initial 
velocities, angular or linear, with respect to any 
axis, Additional control terms, such as flap 


deflection may be included by adding the product of 





a. 
the desired deflection and the appropriate stability 
derivative. Any attitude of the airplane is also per- 
mitted. 

In order to fit the two analogue computers avail- 
able certain assumptions were necessary. The plane was 
assumed to be in horizontal level flight with no initial 
velocities except U, , forward velocity along the 

* waxis. Stability axes are chosen with the undis- 
turbed X-axis alligned with the relative wind and the 
Y eaxis a principal axis. The product of inertia will 
be different from zero. The terms which describe power 
effects are permitted to go to zero by assuming gliding 
flight since these power terms are not essential to the 
dynamic effects which are studied. They merely add 
complexity. 

Solving the resultant equations in terms of the 
highest derivative places them in the form required for 
analogue computer solution: 

(r-1) A= Xn t+ %qg + Xap ur + Xe de barn = rg. — eRe 
(1-8) “rt Yn t+ Wo + Yep 4 aire Bal ceil nut 9e 
(r-a) dr © Zau +£ag + Zu wt + Ee by + 4u- pw 
(xr-i0) = ; Niger aah Pan Ms Nop — Neon ae Me 
aaa 7 Pe Ua) -~qr lve 


< 


ape a ~ ON (IT, 7 a 8 + Poy 
aa - 
(1-12) Bs bemus C Te eS + M. ie 


~pr (IeTe) +rt Ie ->* Le 
Ty 4 ly 
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The assumptions made simplify the equations 
sufficiently to permit use of the analogue computers. 
The inertial coupling terms remain. These products 
are handled by the servo-multipliers and permit the 
airframe studied to move with six degrees of freedom. 

The dimensional stability derivatives are obtained 
from Table V, substituting these values in the equations 
above yields the analogue computer equations for the 
basic airframe: 

Longitudinal Equations 
(aT - 13) we —worslut.,O0L64ur + VR 7 ug. ~ $2,aee 
(m- 14) Urs ~. Out 165 | ~ 2-10 ur +13 By 
+ O10 & a fm py 
(17-15) a* ".000 282 U= Ik, Beg -.034¢ur t 16.8 5, -.00184 ur +. 688 pn T.05nt -.08 
Lateral Equations 
(r-io) We +1 Gp r—. (13V -30Lp +22.4 3, 
“Gon. «Weep = cu + Suze 
(Gr-17) Bt +. 1785p. 0S ay —SS3-p + 6255R + B84 Bp 
+ 0155 RN —-. 85 Gn +. 0185 pa 
(r-18) mR 2 +,01I28V - .190R 7, 0508p ~ 5:46 3 
+1218, t.0323-p 7 ‘Ll8pqa —.032% 4p & 

These equations of motion are entered directly 
in the analogue computers. The dimensional form of 
the stability derivatives permits the use of real time, 
that is, one second of machine time is equal to one 
second of real time. The coefficients, which are the 


dimensional stability derivatives, are entered as 


* . y -- 
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potentiometer settings either directly or, if they 
are greater than one, as a fraction which is sub- 
sequently amplified to full value. The potentiometer 
settings along with the parameters they represent 

are presented in Tables VI and VII. 

The stability derivatives which represent the 
change in a variable due to a control movement are 
based on a one radian control displacement. The 
product of the control movement and the corresponding 
stability derivative represents the forcing function 
which disturbs the airframe from its level flight 
equilibrium. Since the equations are entered directly 
with no machine scaling there is a one to one corres= 
pondence between machine units of 100 volts and unit 
movenent, or in the case of angular control movement, 
unit radian throw. For a 10° control movement a ratio 
is set up which establishes the potentiometer setting 
for a 10° throw. A lesser or greater movement would, 
of course, be represented by a proportionately greater 
or lesser potentiometer setting. An example is calcul- 
ated below to obtain the potentiometer setting represent] 
ing Ms, for a 10° control movement. 

Cea 76.8, Reereseurs THe Erecer 
Cavsep BY & 5S, of | RADIAN. 
\ Macumwe Unir = 100 -volkbs 


ae THE m5 POoTENTIOMETER SETTING 15 


“16.8 
{oo 





VourTs OR .168 
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That is to say that when one machine unit, or 
100 volts, is placed across the M,, potentiometer set 
at .768, the analogue of the change to the pitching 
moment caused by an elevator throw of one radian will 
enter the equations of motion as a forcing function. 
Ten degrees would be this value multiplied by the 
fraction er or e134. The other potentiometer 
settings representing the change caused by control 
movements are similarly computed. For the basic air- 
plane these values appear in Tables VI and VII. 

The computer diagram for the simultaneous solution 
of these six non-linear equations of motion is presented 
as Fig. 9. The inertial coupling terms, which are the 
product or moment of inertia factors and velocity terms 
such aS) qr (Tz -Ty) , are handled by the servo-multi- 
plier units as ateakes in Fags 9. 

A check was made on the analogue computer solutions 
obtained for the basic airplane (Case I) at sea level 
using equations (III=28) of Ref. 1 for the longitudinal 
motion and equations (III-97) for the lateral motion. 
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